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Abstract 

Nuclear collisions exhibit long-range rapidity correlations not present in proton-proton 
collisions. Because the correlation structure is wide in relative pseudorapidity and nar- 
row in relative azimuthal angle, it is known as the ridge. Similar ridge structures are 
observed in correlations of particles associated with a jet trigger (the hard ridge) as well 
as correlations without a trigger (the soft ridge). Earlier we argued that the soft ridge 
arises when particles formed in an early Glasma stage later manifest transverse flow. 
We extend this study to address new soft ridge measurements. We then determine the 
contribution of flow to the hard ridge. 

Keywords: Relativistic Heavy Ions, Event-by-event fluctuations. Two Particle 
Correlations. 



1. Introduction 

Jet correlation measurements at the Brookhaven Relativistic Heavy Ion Collider show 
a striking increase in the yield of associated particles in a narrow range in relative az- 
imuthal angle near (jj fa 1 . This region of enhanced particle production is called a ridge 
because it extends over a broad range in relative pseudorapidity. Interestingly, a similar 
ridge of enhanced production has been reported for two-particle correlations without a 
jet trigger |2|. We seek a common explanation for the jet-triggered hard ridge and the 
untriggered soft ridge based on particle production in an early Glasma stage followed by 
transverse flow [5], [J. 

Correlation measurements of high momentum particles provided the first evidence of 
jet quenching in nuclear collisions [5j [6l [7l |8] . New jet-tagged correlation studies observe 
a high transverse momentum trigger particle and measure the yield of associated low 
Pt particles [TJ [SI UHl IH] • The yield is reported as a function of azimuthal angle (p and 
pseudorapidity 77 relative to the direction of the trigger particle in gold-gold, Au+Au, and 
deuteron-gold, d-|-Au, collisions. A jet peak near rj — (p — likely consists of particles 
from the fragmenting jet trigger, assuming the high pt trigger particle is indeed from 
a jet. The hard ridge appears in Au4-Au collisions as an enhancement of the yield at 
small (j) under the jet peak. The ridge is much wider in 77 than the jet peak, which has a 
gaussian width of ^ 0.2 units, and is absent from d+Au collisions. 

The staggering success of such studies has inspired broader interest in correlations 
of all particle types, with and without triggers. Untriggered correlation studies measure 
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the number of particle pairs over a more comprehensive momentum range dominated by 
typical low momentum particles of < 1 GeV [121 lU ESI [111 IISI ESI ■ Particles in this 
range are primarily thermal and well described by hydrodynamics in Au+Au collisions 
El [3 [S] . The soft ridge is an enhancement of correlations in central collisions that - like 
the hard ridge - is broad in ry = 771 — 772 and centered in the near side, i.e., 77 = = 0. The 
soft ridge is greatly reduced in peripheral Au+Au collisions and absent in pp collisions. 

We have described the soft ridge as a consequence of early-stage correlations in concert 
with late-stage transverse flow |4^. Ours is one of a family of models in which particles 
are initially correlated at the point of production [T71 H EB ESI 1201 lU 122] • Flow then 
boosts the correlated particles into a small opening angle in </>. Flow models have also 
been used to describe qualitative features of the hard ridge [23], [22] . 

The hard ridge has most commonly been described as a consequence of the jet passing 
through the flowing bulk matter produced by the nuclear collision [211 [251 [2S1 [23 [2H1 
[2ni [301 [SD [32] ■ While in principle jets and minijets may contribute to the soft ridge 
[55] , the overwhelming success of hydrodynamics at low pt suggests both effects play a 
role. Of particular interest in this regard is work by Shuryak in which jet quenching and 
transverse flow of the bulk play a combined role |34| . Jet quenching produces a near side 
bias by suppressing the away side jet. Near side jet particles are then correlated with 
other particles from the same jet or transversely flowing bulk particles. 

Important motivation for our work comes from data from the PHOBOS collabora- 
tion, which suggests that the hard ridge and possibly the soft ridge may extend over the 
broad range —4 < 77 < 2 [TTj . Correlations over several rapidity units can only originate 
at the earliest stages of an ion collision when the flrst partons are produced [sL 4^. Hydro- 
dynamics, hadronization, resonance decay, and freeze out can modify these correlations, 
but causality limits such effects to a horizon of roughly from one to two rapidity units. 
Many models referenced above are challenged by this data [35 . 

In this paper we extend the model of Ref. |4] to incorporate jet production and 
address the soft and hard ridges. Long range rapidity correlations and the insight they 
provide on early time dynamics are our driving concerns. We take particle production to 
occur through a Glasma state. Our emphasis is on how computed Glasma correlations 
can affect ridge measurements. In Ref. [T we found excellent agreement with the peak 
amplitude and azimuthal width shown in current Au+Au data. In the next section we 
extend this work here to include Cu+Cu systems. We then extend the model of Ref. |4] 
to address varying pt ranges so that we may address the hard ridge. In Sec. [3] we add a 
contribution of jets following the model of Ref. [34j . We extend that model to compute 
both the strength and azimuthal dependence of the jet contribution to the hard and soft 
ridges. We then combine the flow and jet effects and find that correlations of thermally 
produced pairs constitute a significant contribution to the triggered measurement. We 
then discuss how experiments might distinguish the different contributions. 

2. Glasma Correlations 

The theory of Color Glass Condensate (CGC) predicts an early Glasma stage in a high 
energy collision in which particles are produced by strong longitudinal color fields. As 
nuclei collide, the transverse fields of each nucleus are instantaneously transformed into 
longitudinal fields that are approximately uniform in rapidity. These fields are essentially 
random over transverse distances rt larger than the saturation scale Qj^, where Qs ~ 
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1 — 2 GeV. We can think of such field configurations as consisting of a collection of 
longitudinal flux tubes. Flux tubes are ubiquitous in QCD-based descriptions of high 
energy collisions. In the Glasnia they are closely packed and not strictly distinct due 
to saturation. They are, however, uncorrelated for rt > Ql^^- This is their essential 
feature for this work. The Glasma changes to plasma as particles form from the flelds 
and thermalize. 

In the saturation regime, the number of gluons in a rapidity interval Ay is 

N = {dN/dy)Ay ^ a~^QlR\, (1) 

where Ra is the nuclear radius and as is the strong coupling constant at the saturation 
scale Qs [36]. We understand ([ij as the number of flux tubes K ~ {QsRaY times the 
density of gluons per flux tube (x . The scale of correlations is set by 

_ {N^) - {Nf - {N) 

where the brackets denote an average over collision events. This quantity vanishes for 
uncorrelated gluons, for which multiplicity fluctuations are necessarily Poissonian. In 
Ref. U we argued that the Glasma correlation strength is 7?. cx (K)'^^ — (QgR)^^, and 
found that the Glasma contribution to correlations is 

ndN/dy^as{Qs)-\ (3) 

a result consistent with calculations of Dumitru et al. in Ref. f5]. Equations ([iJ and ^ 
constitute initial conditions for the hydrodynamic evolution of the system. 

These Glasma initial conditions affect the flnal state correlations in several ways. 
First, particles emitted from the same tube share a common origin that is localized to 
a very small transverse area, since Qs ^ Ra- Second, the flux tubes correlate particles 
over a large pseudorapidity range. Some of the flux tubes can stretch across the full 
longitudinal extent of the system at times < 1 fm. These flux tubes rapidly fragment. 
At later times, the particles they produce can be separated by large longitudinal distances 
depending on their momenta. Consequently, subsequent scattering and hydrodynamic 
evolution cannot erase their correlations - they are causally disconnected. Third, the 
strength of the correlation depends on the number of flux tubes. The number of tubes 
depends on the centrality and energy of the collision, as well as the transverse area of 
the tube, all of which, in turn, depend on Qs- Finally, as a result of the common origin, 
particles coming from the same tube must have the same initial radial position and feel 
the same effects from flow, independent if their rapidity. 

During the Glasma phase, flux tubes thermalize into partons and pressure builds 
as the systems moves toward an equilibrated Quark Gluon Plasma. Partons initially 
localized in tubes are now localized in small fluid cells with a uniform azimuthal distri- 
bution of particles. Following a Hubble- like expansion of the system, the transverse fluid 
velocity takes the form 7tVt = Arf. All of partons in a fluid cell are boosted radially 
depending on their initial radial position. Consequentially, partons in any given fluid 
cell gain transverse momentum in the radial direction and the relative angle between any 
two momentum vectors in that cell becomes smaller. Furthermore, fluid cells at a larger 
radial position have a larger flnal transverse velocity and the relative angle between par- 
ton momentum vectors is narrower. This angular narrowing depends only on the initial 
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radial position and the small transverse area of the flux tube source. In our simplistic 
view, all flux tubes are uniform in rapidity and extend to the same longitudinal length. 
This provides for a longitudinally uniform fireball that experiences the same radial flow 
at every longitudinal position. In this way, at every rapidity, angular correlations are 
enhanced in the same way and the initial state spatial correlations are both preserved 
through freeze out and represented in momentum correlations. 

Following [3] we define the momentum space correlation function at freeze out as 



where p2(Pi, P2) = dN/dyid?ptidy2(Ppt2 is the pair distribution, and pi(p) = dN/dycPpt 
is the single particle spectrum. We describe the effect of flow using the familiar blast- 
wave model [371 [Ml [Ml SOI [H]- In this model, the single particle spectrum is pi(p) = 
/ /(x,p) dT, where /(x,p) = (27r)~^ exp{— p'^u^/T} is the Boltzmann phase-space den- 
sity and dT = p^da^ is the differential element of the Cooper-Frye freeze out surface. We 
assume a constant proper time freeze out, so that dT = r^^^mt cosh(y — ri)drjd'^rt, where 
77 — (1/2) ln((t + z)/{t — z)) is the spatial rapidity. 

We argue in Ref. [4 that the flnal-state momentum space correlation function is 



where ni(x) — J /(x,p)(ir. The spatial correlation function c(xi,X2) depends on the 
Glasma conditions as follows 



where rj = r^i — rt2 is the relative transverse position, and Rt = (rti -f- Vt2)/2 is the 
average position. The delta function accounts for the fact that Glasma correlations are 
highly localized to rt < Q^^- The factor p^^CRt) describes the transverse distribution 
of the flux tubes in the collision volume, which we assume follows the thickness function 
of the colliding nuclei. 

We comment that the form of ^ holds as long as 7?. is unmodified from its initial 
Glasma value by particle production and hydrodynamic evolution. This is only strictly 
true as long as a) subsequent evolution doesn't change the relative number of particles in 
the rapidity interval of interest; and b) the number of observed hadrons is proportional 
to the initial number of gluons. Causality prevents these effects from altering TZ for 
truly long range correlations, jryi — 772! > 1 — 2. This would hold for smaller rapidities in 
Glasma theory as long as boost invariance is a reasonable approximation, since dN/dy 
is then a hydrodynamic constant of motion in each event. Moreover, assumption (b) is 
common in Glasma/CGC calculations. On the other hand, for \r]i — t]2\ < 1 — 2, the 
experimental dN/dy is not flat and will change with time due to particle diffusion and 
number changing processes; see |42j . For now, we will assume that TZ is constant and 
defer the hydrodynamic modification for later work. 

The analysis in Ref. [4J focused on a measurement of the near side peak of the soft 
ridge in 200 GeV Au+Au using the observable Ap/^/p = {psib ~ Pre/)/ y/Pref O- The 
quantity psib{4>,vi) represents the distribution of "sibling pairs" from the same event, as 
a function of relative pseudorapidity and relative azimuthal angle, and is comparable to 
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?'(Pl,P2) = /52(Pl,P2) - Pl(Pl)/0l(P2) 



(4) 




c(xi,X2) =TZ5{vt)Pp^{Iit), 



(6) 




Figure 1: Au+Au amplitude as a function of centrality {u = 2Ni,i„/Npart) for both 200 and 62 GeV. 
Solid lines represent Glasma initial collisions TZdN/ dy from (jsjl plus blast wave flow |9]l; parameters are 
unchanged from Ref. 4 . Dashed lines replace Glasma with independent NN initial conditions. 



our p2. The quantity pref{(t>,ri) represents uncorrelated pairs from mixed events and is 
equivalent to the square of our pi. The difference (psib — Pref) is a measure of correlated 
pairs and is comparable to the integral of (|5| over the transverse momenta and average 
azimuthal angle — {(jji + (/'2)/2. It is convenient to compute the quantity 



Pref J Pl{Pn)Pl{Pt2)d'^Ptld'^Pt2 



(7) 



which is independent of the overall scale of the multiplicity. We emphasize that this 
quantity includes correlated pairs in which both particles can have any momentum, 
while measurements of the hard ridge correlate particles from different pt ranges. We 
will extend our approach to address such quantities below. 

To construct the observed quantity Ap/^p^e/, we notice that r(pi,p2) computed 
from eq. ([t]) is proportional to the correlations strength TZ We can write 

-nF{<t>), (8) 



Pref 



where F{(f) is normalized such that J^^ F{(j))d(j) — 1. The distribution F(^) depends 
only on the blast- wave parameters ^m/T and Vs, and represents the angular correlations 
of particles from flux tubes after hydrodynamic expansion. The factor TZ scales the 
strength of the correlations with both energy and centrality and determines the rapidity 
dependence (which is flat in this case). 

We now combine ([t]) and ([T|) to obtain the observed quantity 

^pl^p^f - K-RdN/dy F{<P), (9) 

where we equate the factor TZdN/dy with ([3| , which accounts for all of the Glasma energy 
and centrality dependence. The scale constant k is independent of energy. As described 
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Figure 2: Cu+Cu amplitude as a function of centrality {i/ = 2Ni,i„/Npart) for both 200 and 62 GeV. 
Solid and dashed lines compare Glasma and independent NN initial conditions as in Fig. ^ 



in [3] we set k only for Au+Au 200GeV collisions such that F{(f)) for the most central 
collisions is aligned with the most central data point. Although blast-wave parameters 
have some energy dependence, the Glasma factor ([3| allows for strong agreement with 
the 62 GeV data without further adjustment of k. 

We now apply (|9| and ([3| to address new and forthcoming data for Ap/ ^/Pref and 
the azimuthal widtn Cu+Cu for both 200 and 62 GeV as a function of centrality. No 
parameters are adjusted from [3] - in that sense, these results are predictions. The 
has a dependence on the density of participants Ppart that is determined in Ref. |36j . 
The saturation scale in central Cu+Cu is then 

Q.(cu)^^Q.(Au)^ ^^--;;-'^^-^'^-;; . do) 

Ppart(Au, central) 

After this scaling, we take the relative centrality dependence of Qs to be the same in 
the Cu and Au systems. This assumption can be refined by measuring the centrality 
dependence of dN/dy in Cu+Cu. Similarly, we obtain the blast wave parameters T and 
Vs in Cu+Cu from the Au+Au values by assuming that they scale with the number of 
participants. These assumptions can be refined as single particle spectra measurements 
as a function of centrality become available. 

In Fig. [Tjwe reproduce calculations of Ap/ ^Pref from |1] using (joj) for Au+Au at 
200 GeV and 62 GeV. Preliminary STAR data are from [2]. Here, we compare these 
calculations to dashed curves computed assuming initial conditions from the wounded 
nucleon model rather than Glasma. The wounded nuclcon model implies TZdN/dy in ([9]) 
is independent of centrality. We fix its value to the most peripheral 200 GeV Au data 
as a proxy for pp, and take it to be independent of energy. Note that we might have 
alternatively fixed the wounded nucleon model to peripheral collisions at all energies, but 
this would introduce new unconstrained parameters. The importance of the centrality 
dependence of the Glasma initial condition (|3| is evident, particularly at the highest 
energy. 

Figure [2] shows our prediction for the soft ridge amplitude in Cu+Cu 200 and 62 GeV 
systems as a function of centrality compared to preliminary STAR data ^43) . In both 
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Figure 3: Comparison of the previously published angular width calculations for Au+Au (black line) 
with the angular widths for Cu+Cu systems (blue line, hatched error band) compared with preliminary 
STAR Au+Au 200 and 62 GeV data as well as Cu+Cu 200 GeV data [43|. Width calculations remain 
independent of energy and nearly independent of system 



panels, the dashed Hne represents wounded nucleon model initial conditions as in Fig. 
[T] The dashed lines are included to show how the energy and system-size dependence 
affects the calculation when the Glasma dependence ^ is omitted. The solid lines are 
the result of including the Glasma scaling ^ adjusted by ([To| . The error band represents 
a 10% uncertainty in the blast-wave parameters plus an additional uncertainty in the 
parameterization of Qs that increases with decreasing centrality. We comment that 
the wounded nucleon model with energy independent TZdN/ dy adequately describes the 
data in peripheral collisions but fails in central collisions. This also makes the wounded 
nucleon model prediction larger than the Glasma curve in 62 GeV Cu-|-Cu. 

In Fig. [3] we show the soft ridge azimuthal width in Cu-|-Cu systems compared to 
preliminary STAR measurements as compared to previously published Au-|-Au result 
|4]. We have also included preliminary STAR measurement of the Cu-|-Cu width [43] to 
compare with other preliminary STAR measurements of the soft ridge in Au-|-Au width 
[2]. In [3] we find that the azimuthal width of the near side peak of the soft ridge in 
Au-|-Au is due to radial flow, is constant with a change in energy, and is relatively uniform 
with change in centrality. The error band is representative of the uncertainty of fitting 
an offset gaussian to the angular calculation. Since the azimuthal width is completely 
determined by radial flow, which depends completely on the choice of centrality and 
blast-wave parameterizations, and all of those parameters have remained unchanged, we 
calculate the same enhancement in the width for Cu-|-Cu as Au-|-Au. The black line 
extending to = 6 is the Au-|- Au result from |3j , and the overlaid blue line extending to 
J/ = 4 with the hatched error band is the Cu-f Cu result. Again, as with the Au-f Au result, 
the Cu-|-Cu result is independent of energy since the measured transverse expansion does 
not depend on energy. 

A key feature of the flow-based descriptions of the ridge is that it is the angular width, 
(Jr, of correlated pairs decreases as the mean pt of the pair increases. The greater the 
radial boost given to a fluid cell, the narrower the relative angle between the momentum 
vectors of particles in that cell. A very high momentum correlated pair is more likely to 
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Figure 4: Top panel shows the angular width for Au+Au systems with increasing minimum pt^min limits 
compared with preliminary STAR Au+Au 200 GeV data. The dashed line represents the pt,min=0 
calculation shown in Fig.js] The lower panel shows the azimuthal width for most central collisions vs the 
Pt min limit for both Au+Au and Cu+Cu 200 GeV. Preliminary STAR data for Cu+Cu is also shown 

da- 



have come from a fluid cell that received a very large transverse boost. To study whether 
this effect is present, we compute Ap/ ^Pr&s for pairs of pt > pt,min- As pt,min increases, 
we also expect the amplitude of the ridge to decrease, since it is more difficult to find 
higher pf bulk particles. We therefore compute 



Pref 



Pt2,maj: Ptl.max 

I I '^(Pl,P2) 



Pt2,mirz Ptl,„ 



Ptl.max Pt2,ma,x 

I PliPtl) J Pl{Pt2) 

Ptl.min Pt2,mi,i 



— TZF{(j); Ptl,mimPtl,max, Pt2 ,min , Pt2 ,max) , (H) 

where the integration measures are the same as in ^. To obt ain the measured ratio we 
write 

J'OO 
i , , ^ , 

Ap/Vp^ " Fid,) f^p, ■ ^ ' 



For increasing values oipt^min we calculate the correlation function as before using 
and (12 1 and find that the azimuthal width does indeed decrease as shown in the upper 
panel of Fig. |4] 

The upper panel of Fig. [5] shows the correlation amplitude vs. centrality for different 
choices oi pt^min- The amplitude decreases with increasing pt,min because the number 
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Figure 5: Top panel: Au+Au 200 GeV amplitude calculations for increasing pt,min limits. The dashed 
line is calculation with pt,min = 0. Lower panel: the soft ridge amplitude for most central collisions 
plotted as function of the pt,min limit for both Au+Au and Cu+Cu 200 GeV. Preliminary STAR data 
is also shown |43| . 

of particles contributing to correlations is reduced. The lower panel of Fig. |5] shows 
the amplitude of Ap/ ^Pref for the most central collision as a function of the choice of 
Pt,min for Au+Au and Cu+Cu at 200 GeV. We compare these calculations to preliminary 
Cu+Cu data Similarly, the blue curve in the lower panel of Fig. |4] represents the 
azimuthal width of the soft ridge in most central collisions as a function of pt,min- 

We see that the azimuthal width of the hard ridge is smaller than that of of the 
soft ridge, but as the pt,min limit of the soft ridge is increased, the amplitude of the 
correlations drops and the azimuthal width narrows. In the pt range of the hard ridge, 
it appears that the azimuthal width could be narrow enough, but the amplitude is not 
directly comparable. To understand this difference, we must understand the differences 
in the two measurements. The most significant difference is the choice of the momen- 
tum range of the correlated particles. The hard ridge measurement analyzes the yield of 
associated particles per jet trigger where the associated particle pt range and the trig- 
ger range do not overlap. The soft ridge measurement, however, finds the number of 
correlated pairs per particle where both particles are in the same range with pt above 
minimum bias. The normalization of the soft ridge is found by taking the square root of 
the uncorrelated pair reference spectrum. 

STAR measures the hard ridge, or yield of associated particles per jet trigger, for 
Au+Au 200 GeV for 3 < Pt,trigg < 4 GeV with 2 < pt,assoc < 3 p. Identifying pn 
with the trigger range and the associated range with pt2, we calculate Ap/^p^e/ and 
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transform to yield by 



Yields ^Jffe^l. (13) 



At higher ranges of pti.2 the contribution from jets should become more significant. It 
is important therefore to know the relative contribution of thermal particles and jet 
particles. As will be discussed in more detail later, we decompose the total particle 
spectrum into thermal bulk and jet fractions, and to obtain the contribution of bulk 



correlations to the hard ridge, we multiply (13 1 by the bulk fraction J pi{pti)/ J Ptot{pti) 



where J Ptot{pti) is the total number of particles in the range of pti- 

The blue curve in Fig. [6] represents the contribution to the hard ridge from only 
thermal bulk pairs. As can be seen on the figure, bulk-bulk correlations contribute sig- 
nificantly to the amplitude of the triggered measurement, but seems to have a somewhat 
narrow profile in azimuth. It was shown in |34j that a jet acquires angular correlations 
with fiowing matter due to quenching, but the width of the correlation is wider than the 
data. The contribution of jet correlations with bulk particles could make up the differ- 
ence between the blue curve in Fig. [6] and the data by increasing both the amplitude and 
the width of the calculation. 

In the next section we combine a theory angular correlations from [34 with spatial 
correlations of jets and flux tubes to obtain a jet-bulk contribution to the hard ridge. 



3. Jets, Glasma, and Correlations 

As the Pt of correlated particles is increased, the contributions from jets should be- 
come prevalent, particularly for small 77. At small rapidity differences, correlations of jet 
particles with fragments should be large, but restricted to the size of the jet cone, and the 
transfer of momentum from jet particles to bulk particles is causally limited to ~ 1 — 2 
units in rapidity. The existence of correlations with jet particles at larger 77 would require 
a correlation early in the collision that remains through the longitudinal expansion of 
the system and is still present at freeze out. Both the hard collisions and flux tubes are 
made in the initial moments of the nuclear collision. Assuming that the entire overlap 
region of the colliding nuclei is in the saturation regime, flux tubes would fill the collision 
volume and a jet formed at any transverse position would be accompanied by a flux tube 
at the same position. Since the flux tube extends to large rapidities, the correlation of 
particles from the tube and the jet can extend to large 77. Angular correlations arise 
since particles from the tube acquire a radial trajectory from flow as before, but the jet 
trajectory has a bias in the radial direction due to quenching |Mj. 

We construct a distribution of jets as follows. We assume that jets are produced with a 
hard scattering rate /o(pi) that is independent of position, multiplied by a spatial profile 
Prod oc (1 — ri/R\) that is roughly proportional to the density of binary collisions. The 
phase space density of jet particles is then 

/j(xi,pi) = foiPi)Pprod{ri)S{ri,(j)i), (14) 

where S is the survival probability of a jet due to jet quenching. In practice /o(pi) 
cancels in Ap/pref so we need not specify it. We follow Ref. ^34| and take 

^(ri, 0i) = exp{-L{ri,(l)i)/labs), (15) 
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Figure 6: Angular profile of the jet triggered ridge in a rapidity range away from the jet peak. The 
solid black line combines of long range correlations. Bulk-bulk (blue line) and jet-bulk (dashed red 
line) contributions are shown separately. The jet fraction is determined by ps = 1.25 GeV; bulk-bulk 
correlations make up ~ 75% of the total amplitude. 



where labs — 0.25/m is the jet attenuation length. The survival of the jet depends on 
the path it takes out of the medium 



L(ri, 0i) ^ \jR\-rl siv?{<j)i) ~ n cos(0i). (16) 



The path (16 1 is the distance a jet would have to travel out of a circular transverse area 



at an angle (pi with respect to the radial vector pointing to its position of production ri 



[33]. In view of ( 16 1 and ( 15 ), the shortest path, which is the path a jet is most likely to 
survive, is one that is radially outward from its position of production. Although a jet 
production is at a minimum on the surface, r « Ra, with little material for the jet to 
pass through, the survival probability is maximum in any direction (not pointing into the 
volume) . The resulting angular correlations are weak, since only a small fraction of phase 
space is restricted by quenching. The largest probability of production would occur at 
the center where L « Ra in all directions, but there are no correlations in this case since 
quenching would be maximum. More concisely, jets are less likely to be produced at the 
surface and have a wide angular distribution. The most jets are produced the center and 
would have the narrowest correlation with radially flowing particles, but have the highest 
probability of being quenched. As the production point of the jet moves from the center 
toward the surface, the probability of production decreases, the probability of survival 
increases, and the angular correlation with radially flowing particles widens. Integration 
over all possibilities determines the width of the correlations. 

The calculation follows the analysis of ([S]), but with the first particle from a jet and 
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Figure 7: Comparison of tiie angular shape of tiie jet triggered ridge to long range bulk-bulk (blue line) 
and jet-bulk (dashed red line) correlations arbitrarily normalized to the measured peak height. 



the second from a flux tube, so that 

'^,7b(Pi,P2H/ Cjb(xi,X2 1—^ — -dridT2. (17) 

J "ij(xi) riiB(x2) 

The correlation function c jj^(ii.i,X2) requires that a jet and a bulk particle must come 
from the same radial position. The rest of the equation accounts for the different the 
spectra for jet and bulk particles. We write 

C,JX1,X2)=^^|0^C(X1,X2), (18) 

where c is given by ^ . 

To relate the correlation strength TZj^ to the bulk correlation strength TZ discussed 
earlier, we assume that the hard scattering rate is independent of the flux tube dynamics. 
Recall that the bulk quantity TZ in (2) is related to the number of flux tubes; see ([s]). 
If we take the fraction of jet and bulk particles per flux tube to be independent of 
the number of flux tubes, we can write (Nj) = a{N), where a and /3 may depend on 
momentum, but do not vary event by event. We then follow [ll] to find (N^) — (3{N), 
and {N,N^) = al3{N{N - 1)), so that 

^ {N,N,) ~ {N,){NJ ^ aP{N{N-l))-a{N)p{N) 
{N,){N,} a{N}P{N) 

{NiN^l))-{N)^ 
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Figure 8: Au+Au 200 GeV soft ridge most central amplitude (top panel) and azimuthal width (lower 
panel) vs, the Pt^min predictions with = 1.25 GeV. Blue lines represent the most central points of the 
soft ridge calculation for varying pt^min as shown in Fig, [s] (width) and Fig. [s] (amplitude) . Red lines 
represent similar calculations for jet-bulk calculations, ancTblack lines represent the total as determined 
by the relative bulk and jet fractions. 



Therefore, the addition of jets to the total multiplicity doesn't change the correlation 
strength. In essence, the beam jet associated with the hard process is just another flux 
tube in the high density Glasma state. 

We can now rewrite ([9| for jet-bulk correlations as 

^Pjs/VP^f = ^ndN,et/dy F,M)- (20) 



Our calculation of yield and implementation of lower pt limits follows ( 13 1 and ( |11[ ) but, 
this time, we scale by the jet fraction J pi,j(pti)/ / Pi,tot{pti)- 

In order to compute the amplitude and azimuthal width of the hard ridge in Fig. [6] we 
must determine the relative contributions from both bulk-bulk and jet-bulk correlations. 
The measured pt spectrum follows an exponential behavior at low pt and a power law 
behavior where jets play a larger role, see e.g. Ref. [15]. Our blast wave formulation 
describes the exponential behavior of the low pt spectrum well. The scale Ps at which 
the spectrum begins to deviate from exponential behavior is proportional to Qs in Glasma 
theory, but the proportionality constant is not known. This introduces a free parameter 
- Ps at y/s — 200 GeV - that we flx below. We then find the number of jet particles by 
taking the difference between the total number of particles and the number of thermal 
particles pij = pi^tot — Pi,b- We take pi^tot from the measured spectrum in Ref. [15] 
and pi^B from the blast wave calculation with the appropriate normalization. 

We now calculate the combined effect of Glasma, flow, and jets on the correlation 
function. Adding the bulk-bulk and jet-bulk contributions, we obtain 



(21) 



In order to compare to the yield of associated particles in the hard ridge, we combine 
(21 1 and (13 1 with the appropriate integration limits. We find that the agreement with 
the data in Fig. [6] requires Ps = 1.25 GeV. The dashed red curve in Fig. |6] represents 
the contribution to the yield from correlations of jet and bulk thermal particles. This 
contribution is too wide. On the other hand, the bulk-bulk correlation function describing 
the effect of flow alone, which is given by the blue curve, is too narrow and the computed 
peak height is too small. To emphasize the disagreement of the angular shapes, we show 
the contributions normalized to the peak in Fig. [7] The combination of the two effects 
shown as the black curve in Fig. [6] gives nice agreement with both the amplitude and 
azimuthal width. 

We now compare hard and soft ridge measurements directly by computing the mo- 
mentum dependent correlation function (21 1. Results are shown in Figs. [S] and |9] At 
low pt.min the contribution from jets is negligible, therefore the amplitude and width 
of (21 1 is determined by the bulk-bulk term. As the pt^min is increased, both the am- 
plitude and the azimuthal width of the bulk-bulk term decreases, while the amplitude 
of the jet-bulk term increases. The azimuthal width of jet-bulk correlations is roughly 
independent oi pt,min\ the growth in the figures for pt,min > 1 GeV is due to the growth 
of the jet fraction. Jet-bulk correlations should become a more significant fraction of the 
total as pt,min IS increased, and the azimuthal width of the ridge should increase toward 
the jet-bulk width. The trends shown in Figs. |4] and [9] follow preliminary data from Ref. 
[43] . The small difference may reflect the fact that our calculation omits the jet peak. 

We emphasize that the decrease of the bulk-bulk contribution to ar with increasing 
Pt in Figs. [8] and [9] is a direct consequence of transverse flow and, consequently, is a flrm 
prediction if the jet-bulk contribution is neglected. The role of jets and other phenomena 
like recombination are less clear. We have chosen the model of Ref. 34J because it relies 
only of the well-studied phenomena of jet quenching. Our calculations using this model 
predict that the width would increase for higher pt ranges. 



4. Summary 

We have studied the contribution to the jet-triggered hard ridge and the untriggered 
soft ridge from long range correlations due to particle production from a Glasma state 
followed by transverse flow [3], [4]. To explore the effect of jet production and quenching 
on the ridge we included a model of jet-bulk correlations following [34;. We found that 
transverse bulk flow and jet production affect the transverse momentum dependence of 
the ridge in different ways, as shown in Figs. [8] and [9] If transverse flow is the only 
contribution, we predict that the azimuthal width of the near side peak ar will decrease 
as Pt increases. Jets may introduce new behavior depending on the way in which they 
influence the ridge. Many other effects such as recombination may influence the hard 
ridge. However, we emphasize that we cannot explain the magnitude of the hard ridge 
yield without a substantial soft component. 

We have emphasized the quantitative comparison to data and the information that 
such comparisons can provide. That said, our current model is still schematic and there 
is a lot of data that we have omitted. We have not discussed the rapidity extent of the 
ridge. This is inflnite in our approximate Glasma calculation, but is flnite in experi- 
ments. We will address this rapidity dependence elsewhere, since that analysis requires 
quantum corrections to the glasma as well as viscous corrections to the hydrodynamic 
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Figure 9: Same as Fig. [s] but for Cu+Cu 200 GeV, with data from [iS] . 

treatment [121 146]. In addition, we do not address the away side region, > 7r/2. The 
broad enhancement measured there is attributed to momentum conservation and can is 
described by hydrodynamic calculations [22]. We also have not computed the jet peak, 
which contains information on the interactions of the jet with matter |32|. Further work 
is needed [55] . 

More generally, experiments have treated the soft ridge and the hard ridge as separate 
phenomenon. The soft ridge measurements intend to study properties of the medium, 
while the hard ridge measurement intends to study the effects of jets on the medium. 
The observables used in these studies are different, though not unrelated. We feel that 
the experimental situation may be clarified by using common techniques and observables 
in the different regimes. 

The fundamental reason for studying these complex phenomena lies in the fact that 
long range correlations provide information on the very earliest stage of particle produc- 
tion. The overall amplitude of correlations is fixed by the CGC-Glasma dependence of 
TZdN/dy. We have seen that this behavior explains the energy dependence and system 
size dependence of the soft ridge. Moreover, as the pt range of correlations is increased 
toward jet-dominated part of spectrum, long range correlations continue to influence the 
near side peak because flux tubes induce correlations of jets with low pt bulk particles. 
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